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Abstract. We consider the velocity fluctuations of a system of particles described 
by the Inelastic Maxwell Model. The present work extends the methods, previously 
employed to obtain the one-particle velocity distribution function, to the study of the 
two particle correlations. Results regarding both the homogeneous cooling process and 
the steady state driven regime are presented. In particular we obtain the form of the 
pair correlation function in the scaling region of the homogeneous cooling process and 
show that some of its moments diverge. This fact has repercussions on the behavior 
of the energy fluctuations of the model. 



Velocity fluctuations in a one dimensional Inelastic Maxwell model. 



2 



1. Introduction 

In recent years the understanding of the physics of granular materials has taken great 
strides. The dynamical properties of particles experiencing mutual inelastic collisions 
has been thoroughly studied experimentally, theoretically and by computer simulation. 
Such an effort has lead to the discovery and to the formulation of new phenomena and 
properties pp. Among these properties a special place is occupied by the homogeneous 
cooling state (HCS), i.e. the state achieved by a granular gas, initially in motion, under 
the effect of the energy loss caused by inelastic collisions. Loosely speaking the HCS 
plays for granular gases a role analogous to the Maxwellian for molecular elastic gases. 
Although the properties of the HCS are known in detail, its explicit form can be obtained 
as series expansion only for some specific models such as the inelastic hard-sphere model 
(IHS) 0. 

The prototype model for the study of granular systems is represented by an 
assembly of smooth inelastic hard spheres, characterized by a constant coefficient of 
normal restitution. For such a model various authors have derived the Boltzmann and 
the Boltzmann-Enskog equations describing the evolution of the reduced one-particle 
velocity distribution [3]. However, since these approaches remain mathematically hard 
to solve, a simpler mathematical model, the Inelastic Maxwell Model (IMM), where 
the collision rate between the particles is assumed to be independent of the relative 
velocity of the colliding pair, has been put forward. In this model the spatial structure 
is neglected and only the velocity of the particles specifies the state of the system. The 
IMM is nevertheless useful and studied because it lends itself to analytical solution 
in one dimension, thus providing a benchmark to test approximate treatments [I]. 
In the homogeneous free cooling case [5], the evolution equation for the velocity 
distribution has a scaling solution that can be expressed in an analytical closed form, 
with high energy tails described by an algebraic decay: the exponent does not depend 
on the restitution coefficient. Moments of the velocity distribution exhibit multiscaling 
asymptotic behavior [6]. The Inelastic Maxwell Model is quite simplified with respect to 
inelastic hard spheres and other realistic models of dilute granular materials, nevertheless 
in the past it has been considered an important starting point for granular kinetic 
theories [7j. As stated by Ernst and Brito [8]: "What harmonic oscillators are for 
quantum mechanics, and dumb-bells for polymer physics, is what elastic and inelastic 
Maxwell models are for kinetic theory." 

Most of the literature on the kinetic theory of granular gases focuses on the 
single particle distribution function. This is in analogy with the relevance that the 
Molecular Chaos approximation has for molecular (i.e. elastic) gases. On the other 
hand, the inelasticity of collisions in granular gases makes this assumption more delicate: 
numerical and experimental evidences show a stronger tendency of granular systems 
to enhance correlations, often appearing in the form of spatial structures [9] [TO], [TTj . 
Fluctuations have been investigated by various authors P21 [13l [T4j . However, only 
recently the two-particles distribution function has come under scrutiny, in particular by 
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Brey and coworkers who considered its application to the study of the energy fluctuation 
in the homogeneous cooling state of inelastic hard spheres [15] . Their study focuses on 
the effect of inelasticity on the 1/N deviations from Molecular Chaos. Here, our aim 
is to apply similar analysis to the one dimensional Inelastic Maxwell Model. This is 
interesting because, with a few controlled approximations, one obtains the asymptotic 
pair correlation function in a closed form, and all time dependencies of its two-particles 
velocity moments, getting further than the original work of Brey et al., where only the 
asymptotic moments, in particular those required to calculate energy fluctuations, were 
explicitly obtained. 

This paper is organized as follows. In section [2] the evolution equation for the IMM 
is presented and the equations for the various distribution functions introduced. In 
section [3] the dynamical equations are solved for the moments of the single and two- 
particles distribution functions. The asymptotic scaling state is discussed, for the pair 
correlation function, in section SJ. Finally, in Section [5] the effects of an external driving 
is considered and in [6] the concluding remarks are presented. 



2. Evolution equations for the distribution functions 

We consider a system of N particles, each characterized by a scalar velocity Vi, with 
...,N. The Inelastic Maxwell Model assumes that the state T = (vi,v 2 , ..,Ujv) is 
modified by elementary collision events, realized by changing the velocities (vi, vj) of a 
randomly selected pair of particles according to the rule: 

v'i = JVi + (1 - j)Vj (1) 
Vj = (1 - j)Vi + JVj 

where 7 = and a is the coefficient of restitution. The system cools down because 
in each collision an amount AE of kinetic energy, given by 

AE = -^(l-a*)(v>-v>y, (2) 

is dissipated, where m is the mass of a particle. Since the IMM is not endowed with a 
spatial structure such a cooling process is homogeneous. An observable A(T(t)) evolves 
according to 

A(T(t)) = exp(t£A(T(0))) (3) 
where the generator C is 

1 N N 

1=1 j^i 

and the operator T acts on an arbitrary function S(vi, Vj) of the velocities of particles 
% and j in the following way: 

T(v h Vj)S(v h Vj) = Siv'i, v'j) - S(v h Vj) (5) 

The "time" t is a collision counter and represents the clock of the model. 
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Following closely the derivation presented by Brey et al.[T5|, in order to set up the 
evolution equations for the system, we introduce the following distribution functions 

N 

F 1 (u 1 ,t)='£(6(u 1 -v i {t))) (6) 

i=l 

N N 

f 2 ( Ui , u 2 , t) = Yl Y,^ ui - u i(*))%2 - v j (*))> ( 7 ) 

i=l j^i 

N N N 

F 3 (u u u 2 , u 3 , t) = ( 5 ( Ul - < t )) 6 ^ - Vj(t))6(u 3 - v k {t))) (8) 

i=l jjki k^i,j 

where (■) stands for an average over an ensemble of trajectories with different initial 
conditions (in section [5j where the effect of a thermal bath will be considered, this will 
be an average over realizations of the noise). 

A hierarchy of equations can be derived for these distribution functions, whose first 
two equations read: 

d f 

Fi («i,t)= / du 2 T(u 1 ,U2)F 2 {u 1 ,u 2 ,t) (9) 



dt 
and 

d 



F 2 (wi, « 2 , *) = T(ui, u 2 )F 2 (u 1 , u 2 , t) + 
dt 

+ J du 3 [T(ux,u 3 ) +T(u 2 ,u 3 )]F 3 (u 1 ,u 2 ,u 3 ,t) (10) 

where the inverse binary collision operator, T, is defined for a generic function S(ui,Uj) 
of the velocities by the rule: 

T(ui,u j )S , (n i ,n,-) = — S(u*,u*) - S(u h Uj) (11) 
a 

which transforms the velocities [u^Uj) into their pre-collisional values u*,u*, obtained 
by inverting eq. ([T]). Following a standard statistical procedure we consider the following 
decompositions of the distribution functions: 

F 2 (« 1 ,Ti 2 ,t) = F 1 (« 1 ,t)F 1 (u 2 ,t) + G 2 («i,« 2 ,t) (12) 

F 3 (u!,u 2 ,u 3 ,t) = F 1 (ui,t)F 1 (u 2 ,t)F 1 (u 3 ,t) + G 2 (iii, u 2 , t)Fi(u 3 , t) 
+ G 2 (u 1 ,u 3 ,t)Fi(u 2 ,t) + G 2 (u 2 ,u 3 ,t)F 1 (u 1 ,t) 
+ G 3 (ui,u 2 ,u 3 ,t) (13) 

After substituting these expressions into and ffTUl) and dropping the term containing 
G 3 we obtain a pair of closed equations for F 1 and G 2 . Before proceeding further we 
also define, for later convenience, the following normalized distributions: 

h(u lt t) = ^F 1 (u li t) (14) 
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h 2 (u u u 2 , t) = — G 2 {u x , u 2 , t) (15) 
which obey the sum rules 

dmhfat) = 1 (16) 

duidu 2 h 2 (ui,u 2 ,t) = —1 (17) 

du 2 h 2 (ui,u 2 ,t) = -/i(iti,i) (18) 
Hence eqs. (jSJ) and ffTUl) can be rewritten as 

^:/i(mi,t)= / du 2 T{u 1 ,u 2 )[f 1 {u 1 ,T)f 1 {u 2 ,r) + — h 2 (ui,u 2 ,T)} (19) 




and 



—h 2 (ui,u 2 ,r) = r(ui,M 2 )[/i(Mi,*)/i(u2,r) + -^/i 2 Cui,u 2 ,t)] 
ar iv 



+ y dttsr(ui,U 3 )[/l2(ui,U2,r)/i(M3,r)+/l2(M2,M3,*)/l( u l> r )] 

+ / du 3 T(ii 2 ,M 3 )] [^2(^1,^2, r)/i(w 3 ,r) +/i 2 (Wl,«3,T)/l(w 2 ,T)] (20) 



where we have redefined the time variable r = iVt. 

In order to solve (fT9l) and fl20|) we slightly generalize the method, originally 
introduced by Bobylev [16] and consider the following Fourier transforms of the 
distributions fi and h 2 



h{h,t) = J d Ul e lk ^ h{u u t) (21) 

h 2 {h,k 2 ,t) = J J ' du x du 2 e ik ^ +ik ^h 2 {u x ,u 2 ,t) (22) 

The function h 2 (ki,k 2 ,t) is symmetric and has the property h 2 (ki,0,t) = — fi(ki,t) as 

a consequence of the sum rule ffTSl) . Substituting these expressions into fTT9~l) and ([2"D1 
we find: 

A/i (fci, r) = A(7*i, r)/i((l - 7 )*i, r) - A(*i, r)A(0, r) 

+ -^(7*1, (1 - 7)^1, r) - h 2 (h, 0, r)] (23) 



and 



—h 2 {k 1 , k 2 ,r) = ^-[^2(7^1 + (1 - 7)^2,7^2 + (1 - 7)^1,^) - ^2(^1, h,r)] 
dr N 

+ A ( 7 *! + (1 - 7 )fc 2j r)A( 7 A; 2 + (1 - 7)^1, r) - A(*i, r)A(* 2 , r) 
+ (1 + P 12 )[/i((l - 7)^1, ^(7*1, k 2 , r) - A(0, r)^ 2 (A:i, k 2 , r) 
+ A (7*1, ^(fe, (1 - 7)^1, r) - A(fci, ^Mfe, 0, r)] 

(24) 
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where P\ 2 exchanges the index 1 and 2. 

It is possible to connect some elements of the pair distribution function to observable 
properties. To this purpose we consider the distribution functions, Ud(V) and H a (W), 
of the difference of the velocities V = (ui — u 2 ), and of the sum W = {u\ + u 2 ), 
which are obtained by marginalizing the distribution f 2 (u±, u 2 ,t) = fi(ui,t)fi(u 2 ,t) + 
h 2 (u 1 ,u 2 ,t)/N according to the transformations: 

MV) = j j d Ul du 2 f 2 (u u u 2 )5(V - ( Ul - u 2 )) (25) 

and 

K.(W)= [ [ du 1 du 2 f 2 (u 1 ,u 2 )6(W -(ui + u 2 )) (26) 



We take, now, Fourier-Bobylev transforms of both distribution functions and find 

fl d (k, r) = h(k, r)A(-A;, r) + ^h 2 (k, -k, r) (27) 

Similarly 

U s (k, r) = h(k, T)h(k, t) + ^h 2 (k, k, t) (28) 

Of course, the correction is of order 1/N and vanishes for infinite systems. 

In the following we shall assume 1/N << 1 and drop the corresponding terms in 
(123|) and (|24j) (see Appendix for a discussion of this approximation). Hence, eq. (|23l) 
reduces to the standard equation of the one-dimensional IMM [B] and decouples from 
the evolution equation for h 2 . 

3. Power series solution. 

The distribution functions can be expanded into their moments as follows: 

M n (r) = J dui</i(ui,r) (29) 

Qmn(r) = J J duidu 2 u™v%h 2 (u x ,u 2 ,T), (30) 

obtaining: 

/i^rl^^fM^r) (31) 

n=0 

h 2 (k u k 2 ,r)= jr ^f ^f QUr). 

m,n=0 

(32) 

Inserting these expansions into eqs. (1231) and fpH|) we first recover the moments, Mj, 
evaluated by Ben-Nairn and Krapivski [6] and given by: 

M (r) = 1 (33) 
M 1 (r) = (34) 
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M 2 (r) = M 2 (0)e~ a2T (35) 

M 3 (r) = M 3 (0)e- a3T (36) 

M 4 (r) = [M 4 (0) + 3M 2 2 (0)]e^ a4T - 3M 2 2 (r) (37) 

with the coefficients given by a n = 1 — (1 — 7)" — 7™, a 24 = 67 2 (1 — 7) 2 and £ = 7(1 — 7). 
Notice also that a 2 = 2£ and a 3 = 3£, and a 4 — 2a 2 = — 2£ 2 . In addition, we obtain 
the moments of h 2 using the conditions that the initial velocities are indipendently 
distributed (with a constraint on the total momentum M\ (r) = 0) and the one-particle 
distribution is even: 

Qno(r) = Qon(r) = -M n (r) (38) 

Qni(r) = Qi„(t) = -M n+1 (r) (39) 

Q 22 (r) = [Q 22 (0) - M 4 (0)]e- 2a2T + M 4 (r) (40) 

We can, now, compute the energy fluctuations since 



Nm 



(E 2 (t)) - (E(t)) 2 = ^~{J J duidu 2 u\ulh 2 {u ll U2,T) + J du^f^ii!, r)| 

■{g 22 (r) + M 4 (r)}, (41) 



4 

having defined the total energy as E(t) = y X/ili ^l 7 ")- Recalling the kinetic definition 
of granular temperature T g (r) = mM 2 (r) = 2{E{r)) /N, one has that 

Mll^^^^e^CV), (42) 

with A = Q22 ^2(q* 4 ^ — 6 and B = 2 M ^°^2^ 2 ^ ■ Therefore the energy fluctuations 
decay at a slower rate than the square of the average energy. The situation is analogous 
to what happens to the fourth moment of the distribution function, which also diverges 
if rescaled by the square of the second moment. On the other hand, we notice that the 
energy fluctuations, scale proportionally to N, as in non critical systems: this means 
that a thermal capacity can always be defined, but it grows with time. This is different 
from what happens in the homogeneous cooling of inelastic hard spheres, as discussed 
by Brey et al.[15j. where the ratio between energy fluctuations and the square of average 
energy is constant in the HCS scaling state. 



4. Fluctuations around the scaling solution. 

It is well known that eq. ( 123]) for large N, possesses a scaling solutionjl], where the only 
time dependence occurs via the combination q\{r) = kiv (r), i.e. fi(ki, r) = £o(?i)> with 
^ ; o( r ) = a/ M 2 (t) the thermal velocity. Using such a variable the evolution equation for 
the distribution function takes the scaling form 

- Cgi-^ofe) = 6)(7<7i)£o((l - 7)31) - UnMo) (43) 
aqi 

which has the solution 

£o(gi) = (l + ki|)e- 1911 . (44) 
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Its small q\ singularity reflects the inverse power law tails of the corresponding velocity 
distribution function [8], 4>o(c) which is obtained by applying the inverse Boby lev- Fourier 
transform to eq. ()44|) with the result: 

Mc) = ~ h ~ 212 , (45) 

7T [1 + C z \ z 

with c = u/vo(t). Note that the complete time-dependent velocity distribution reads 
fi(u\,t) = 4> [ui/v (t)]/v (t). We wish to consider, now, the fluctuations around the 
scaling solution. We first define the functions 6 with i = 1, 2 defined as 

6(9i) = (\qx\ + 9i> _kl1 , 6(<?i) = (46) 
and the linearized Maxwell-Boltzmann operator Ai = A(qi; £0(91)) as: 

Ai^(go = c?i^-^(?i) + e ((i - 7)91)^2(791) - eo(o)v(?i) 

+ £o(79iM(l - 7)9i) - 6(51)^(0). (47) 

One can see that 

Ax£ = C6, Ai6 = C6, Ax6 = (48) 

and conclude that £1 and £2 are eigenfunctions of Ai corresponding to the eigenvalues 
C and respectively, whereas £0 is not eigenfunction of Ai. Interestingly, the Bobylev- 
Fourier transforms of £2 and £1 read, respectively: 

02(c) = ^( C( M C )) 

^ = l\ih +Uc) - (49) 

Remarkably, these two eigenfunctions have a similar structure to that of two of 
the three eigenfunctions found by Brey and co-workers, see Equations (65) in their 
paper [15]. This similarity is however incomplete, since in [15] the eigenfunctions of the 
linearized Boltzmann operator were identified to be the hydrodynamic modes. Up to 
our knowledge, a study of hydrodynamic spectrum for the Inelastic Maxwell Model is 
missing (it has been performed, for Inelastic Hard Spheres in [23]) and therefore we are 
not able to make a similar connection, neither to find the third eigenfunction necessary 
for completing the analogy. 

In order to determine the pair correlation function we, now, rewrite eq. (|24[) in the 
scaling form 

[Ai + A 2 ]x 2 (gi,fc) =U(q 1) g 3 ) (50) 
with X2O71, 92) = h(h, k 2 , r), 

U(qi,q 2 ) = -£o(79i + (1 -7)92)£o((l -7)9i +792) +£o(9i)£o(92) (51) 

and having defined the operator A2 as identical to Ai but acting upon the variable q 2 . 
If 9i92 > formula (15T1) can be cast in the form 

U(qx, 92) = - C[6(9i)£o(9 2 ) + £0(91)6(92) (52) 
+ £1(91)6(92) +6(91)6(92) - 26(9i)6(92)] 
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so that we easily find a solution 

X2(<?i,<?2) = %i 92) (-6(91)6(92) +Ci(9i)Ci(?2) +£226(51)6(92) 

- [6(91)6(92) +6(92)6(91)]) (53) 

where 9(x) is the Heaviside step function and C 22 is an arbitrary constant which may 
be fixed using the boundary conditions. For gig 2 < we could not find an exact 
solution. The most reasonable path to reach an expression for ^2(91; 92) in this part of 
the qi, q2 plane is therefore to expand the r.h.s. of Eq. ( 1501) on the basis 6 (see Appendix 



Appendix A for details) and assume 

2 

Xlfal, 92) = #(-9l92) ^2 C mn€m(Ql)€n(Q2)- (54) 

m,n 

This procedure, which is justified only as an approximate "continuation" of the solution 
in the q±q 2 > region, and will be checked for consistency at the end of this section, 
leads to the following equation for the coefficients C rnn : 

2 2 

[Ai + A 2 ] ^2 Cmn6n(9i)6i(92) = ^ r mn 6ra (91)6(92) (55) 

m,n m.n 

whose approximate solution can be found by expanding T mn in powers of £ up to order 
( 3 (i.e. for small inelasticity). The final result reads: 

O)o = -1 + ^C + 45C 2 (56) 

C i = C 10 = -^C-33C 2 (57) 

di = -l + yC + ^C 2 (58) 
61 

C12 = C 21 = 1 - —C - 128C 2 (59) 

C 02 = C 20 = y C + ^C 2 (60) 

Also in this case C22 results arbitrary. Notice that the form of the pair correlation 
function ^2(91, 92) in the region qiq 2 < now depends on the inelasticity through £ and 
therefore is not universal as the single-particle distribution. 

In the previous section we have seen that the rescaled energy fluctuations, related 
to the moment Q22, diverge as r — > 00. The situation is similar to that encountered in 
the study of the single particle distribution function, where the exponential increase of 
the rescaled fourth moment M4(t)/M 2 (t) 2 was the signature of the fact that the fourth 
moment of 0o( c ) diverges. Is the behavior of Q22 the fingerprint of a similar behavior 
of X2(qi, 92)? Indeed, the size of the energy fluctuations is controlled by the small 
(91, 92) singularities of ^2(91, 92), because these determine the high velocities tails of the 
correlations 0[5]. In order to expose the presence of the tails of X2, we isolate the most 
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singular contribution to it, namely the term proportional to £i(<7i)£i(<72) i n ec L s - dSSD 
and fl54j) 

xT 9 (Qi,<l2) = [Oiqm) + ^(-gig2)C , n]^i(gi)G(g 2 ) 

^ [0(giga) + ^(-gig2)C , 11 ]|g 1 ||g 2 |e- |9l| - |92| 7 (61) 

After Fourier transforming to velocity space we realize that the pair correlation function 
for large value of its arguments decays as: 

1 



FT[ X S 2 



singi 



2^2 



CTC 



(62) 



1^2 



so that the moment c\c\ diverges. Such a result is the counterpart of the divergence of 
the fourth moment of the single particle distribution function. 

Before closing the present section we wish to comment the fact that the projection 
introduces an error due to the truncation of the expansion of U(qi,q 2 ) (see eq.f lA.7l) ). 
How reliable is such an approximation? In order to check the error we computed, for 
various values of a, the following quantity: 

r ~ 2 1 2 

J Jdqidq 2 U(q 1 ,q 2 ) - Em n =o T mn^m{qi)^n(q2) 

L (63) 



/ / dq 1 dq 2 U(qi, q 2 ) 

which represents a measure of the relative error. As we see from fig. [T]the approximation 
becomes poorer and poorer as a — > 0. However, for not too low inelasticities the 
approximation is reasonable and in our opinion this justifies the above procedure, in 
particular Eq. (154")) . for value of a > 0.5. 




Figure 1. The error A of cq. ([63)) as a function of a 
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5. Driven system 

Now, let us consider a system driven by an external Langevin heat-bath which has been 
considered by several authors [T71 HH [191 120] • The velocities of the particles, evolve 
between collisions according to an Ornstein-Uhlenbeck process: 

^ = -r«,(f) + 6(f) (64) 

with 

(&(*)> = 0, mMt')) = ZDSijSit - f) (65) 
The resulting equations for the distribution functions are: 

±f 1 (k 1 ,r)= -(Dkl + Fh^hiKr) 

+ -{hinK r)/i((l - 7 )*i, r) - A(/ci, r)A(0, r) 

+ ^Mlh, (1 - 7)*!, r) - h 2 (h, 0, r)]} (66) 

where an arbitrary mean free r c time has now been introduced for dimensional reasons, 
and 

^-h 2 (k u k 2 , t)= - [D(kf + k\) + r{k 1: £- + k 2 4r)Mh, k 2 r) 
dr dk\ dk 2 

+ -{4fe(7^i + (1 - i)h%, jh + (1 - 7)Ai, r) - ^(/d, fc 2 , r)] 

+ / 1 ( 7 fc 1 + (1 - 7 )fc 2 , r)/ 1 ( 7 fc 2 + (1 - t)*i, r) - A fa, r) Afa, r) 
+ (1 + Pia) [A((l - 7)^1, r)/i 2 (7fc 1 , fc 2 , r) - A(0, r)h 2 {k u k 2 , r) 
+ fihh, r)h 2 (k 2 , (1 - 7 fa, r) - Afa, r)ft 2 fa, 0, r)]} 

(67) 

We look for steady state solutions, setting the time derivatives to zero. By slightly 
modifying the method employed to derive the moments in the cooling case we obtain 
the moments in the non equilibrium steady state regime. We find: 



M = 1 (68) 
M 2 = ^ (69) 

rv c + 7(i- 7 ) 1 ; 

_ 12L>r c M 2 + 6 7 2 (l-7) 2 M 2 2 
4 4rv c + 1 -7 4 - (1 - 7 ) 4 
= 1^ + 7(1-7)1 + 6^(1--,)- 

4rr c + l-7 4 -(l-7) 4 2 V ; 

The granular temperature T g = mM 2 is obtained from equation (16"91 . yielding in the 
elastic case (7 = 0) T g = mD/T as expected. Furthermore, the last equations in the 
elastic case becomes M 4 = 3M|. 
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Similarly we find the coefficients of the moments of the pair correlation h 2 using in 
eq. (167|) the expansion fl30l) 

QnO = Qon = ~M n (71) 

7 ^ ~~ 7 ^-M 2 (72) 



'ii 



Q12 = Q21 = (73) 
7(1 - 7)[1 ~ 27(1 ~ 7)](3M| + M 4 ) + 3T7egfi + 6Jt c <^ 

Ql3 " Qsi " 4rr c + 3 7 (1 - 7) ^ 

7 2 (1 - 7) 2 (3M| + M 4 ) - 2Dt c M 2 



Q 



22 



2IY C + 2 7 (1 - 7) 

_ W1 - ^ 

M 2 2 , (75) 



7 2 (l-7) 2 (C + 3)-2rr c -2 7 (l-7)„ 2 



2Tt c + 2 7 (1 - 7) 

where in the last passage we have introduced the constant C = M4/M. 2 . Again one can 
wonder the ratio between energy fluctuations and the square of granular temperature, 
obtaining: 

(E 2 (t)) - (E(t)) 2 = jV 2rr c + 7(l-7)[2 + 7(7 + 5)] 
T 9 (t) 2 2 2rr B + 7 (l-7)[2-7(l -7)] 

which yields the value N/ 2 in the elastic case [22] . 

Switching to the reduced variable c 2 = v 2 /(2M 2 ) we can look for an expression of 
the distribution function in terms of Sonine polynomials: 

A(c)^e- c2 [l + S2 S 2 (c 2 ) + ..] (77) 



with s 2 = — 1 + M 4 /3M| and S 2 (c 2 ) = | — |c 2 + |c 4 . In practice, one approximates the 
series f J77|) with a finite number of terms and since the leading term is the Maxwellian, 
the closer the system to the elastic limit, the less term suffice to describe the state. 
In the same spirit we assume the following expansion for the two-particle distribution 
function expression 

h 2 (ci, c 2 ) = - /i(ci)/i(c 2 ) + - exp \-{c\ + c 2 )l [A^c-2 

7T 

+ A 2 {\ - c{){\ - c\) + A 3 ( Cl cl + c\ C2 )] (78) 
where the coefficients Ai satisfy the relations: 

= ^ " ^ (79) 



M 2 M, 



A >=- 2 -k+ 2 M (81) 

A straightforward computation shows that Ai — > (for i = 1,2, 3) in the elastic limit 
a — > 1 and in the Brownian limit Tr c — > 00, i.e. when the collision rate is so small that 
grains thermalize with the external bath. 
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6. Conclusion 

We have shown that if the number, N, of particles, experiencing inelastic collisions 
described by the Inelastic Maxwell Model, is finite it is possible to observe correlations of 
order 1/N among the velocities of different particles. Such correlations have been studied 
in two relevant situations: the homogeneous cooling state and the steady state obtained 
by applying a stochastic driving to the system. In the first case we have obtained the 
velocity correlations by solving to order 1/N the equations for the moments of the one 
and two-particles distribution functions which show that the energy fluctuations decrease 
slower than the squared energy. In addition, we have studied the velocity pair correlation 
function in the scaling regime where the one-particle probability distribution is given by 
fx(u,t) = 2[7rv (t)]~ 1 [l+(u/vo(t)) 2 }~ 2 . For small inelasticity we have obtained its explicit 
expression. Interestingly, such a solution shows that the moment Q22 of the velocity 
pair distribution function diverges. We may conjecture that such tails, which are the 
fingerprint of the Maxwell model, will persist in the many-particles correlation functions 
of higher order. These could be in principle computed using the same methods discussed 
above, although the effort required to carry out the program could be exceedingly heavy. 

Finally, we have obtained, by the series expansion method, the pair distribution 
function when the system is subjected to a Langevin driving. In this case the moments 
of the pair correlation are finite up to the fourth order and we believe that the higher 
moments will also be finite. 

As final remark we would like to comment that although the Maxwell model is 
somehow artificial and does not describe any real granular material it offers, as our paper 
illustrates, the possibilty of exploring new aspects of non equilibrium statistical systems. 
Acknowledgments. - U.M.B.M. acknowledges the support of the Project COFIN-MIUR 
2005, 2005027808. 

Appendix A. Projection technique 

Since we are not able to find a solution in the full space we resort to an approximate 
method in the remaining Fourier space. The method consists of projecting the term 
onto the subspace spanned by the function £ n . For the sake of simplicity we define the 
scalar product between two functions / and g 




(A.l) 



and introduce an orthogonal basis, R n (q), of the form 



n 




(A.2) 



1=1 



The orthonormalization conditions for qiq 2 < give the following relations: 




(A.3) 
(A.4) 
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i? 2 (|g|) = 2 1 / 4 (e -56 + 76) (A.5) 
in compact form: 

Rm(\<l\) = J2 Mm ^n (A.6) 

n 

We define, using the Heaviside function 9(x), ^(91,92) = [1 _ 9(qiq 2 )]U(qi, 92) and 
expand with respect to the RnS 

2 

U(qi, to) = Yl U abRa(qi)Rb(q2) + K(q u q 2 ) (A.7) 

a,b=0 

where 

/+00 p+00 
/ dq 1 dq 2 R a {q 1 )R b {q 2 )U(q 1 ,q 2 ) (A.8) 
-00 1/ —00 

and K(qi,q 2 ) represents the part of the function U(qi,q 2 ) orthogonal to the subspace 
spanned by the three functions above. 

Inserting the ansatz ^2(91, 92) = J2 a ,b ^066(91)6(92) into (El), neglecting the term 
K(q 1 ,q 2 ) and using fl A.6|) we obtain 

2 

[A X + A 2 ] ^ CrtUqiMto) = T mnU(Ql)UQ2) (A.9) 

a,b m,n=0 

where T mn = Ylab=oUabM am Mb n - Substituting the expansion of T mn up to third order 
in £ we find the following equation for the coefficients C mn : 

C [C016 (91)6(92) + C o6(9i)6(9 2 ) + C106 (91)6(92) + 2C u 6(9i)6(92) 
+ (C10 + C 12 )6(9i)6(9 2 ) + C 00 6(9i)6(9 2 ) + [On + C 21 ]6(9i)6(9 2 ) 
+ [C02 + C 20 ]6(9i)6(9 2 )] = 7C 3 6(9i)6(9 2 ) 

- (^C 2 + 33C 3 )(6(9i)6(9 2 ) + 6(91)6(92)) + (-2C + 21C 2 + 135C 3 )6(9i)6(9 2 ) 

+ (-C + ^C 2 + 45C 3 )(6(9i)6(9 2 ) + 6(91)6(92)) 
+ (C - 32C 2 - 161C 3 )(6(9i)6(9 2 ) + 6(92)6(91)) 

+ (39C 2 + 147C 3 )6(9i)6(9 2 ) (A.10) 
whose solution is given by eqs. (!56l) - (l60l) 

Appendix B. Small terms 

In order to validate our assumption of neglecting the 1/N terms, we have evaluated 
the terms on the r.h.s. of eq. ( 123]) using the relations given by eqs. ( I37|) -(l40l). We 
obtain that the 1/N terms are null up to third order and the corresponding k\ terms 
are explicitly 

« ... + 1{C 2 [M 4 (0) + 3M 2 2 (0)]e- a ^ - 2CM 4 (r)}fc 1 + ... (B.l) 
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— hilh, (! - t) - h 2 (ki, 0, r) 

« ^{C 2 [Q 22 (0) - M 4 (0)]e- 2a - + 2C 2 M 4 (r)}^ + ... (B.2) 

Comparing both coefficients we can state that our assumption is accurate at least up to 
fourth order. 

[I] Poschel T and Luding S (eds.), Granular Gases, 2001 Lecture Notes in Physics, Vol. 564, Springer, 

Berlin 

[2] Brcy J J, Dufty J W and Santos A, Dissipative dynamics for hard spheres, 1997 J. Stat. Phys 87 
1051 

[3] van Noijc T P C and Ernst M H, Velocity Distributions in Homogeneously Cooling and Heated 

Granular Fluids, 1998 Granular Matter 1 57 
[4] Baldassarri A, Marconi U M B, and Puglisi A, Influence of correlations on the velocity statistics 

of scalar granular gases, 2002 Europhys. Lett. 58 14 and Kinetic Models of Inelastic Gases, 2002 

Math. Mod. Meth. Appl. Sci. 12 965 
[5] Marconi U M B and Puglisi A, Mean-field model of free- cooling inelastic mixtures, 2002 Phys. Rev. 

E 65 051305 

[6] Ben-Nairn E and Krapivsky P L, Multiscaling in inelastic collisions, 2000 Phys. Rev. E 61 R5 
[7] Santos A, Transport coefficients of d-dimensional inelastic Maxwell models, 2002 Physica A 321 
442 

[8] Ernst M H and Brito R, High- energy tails for inelastic Maxwell models, 2002 Eurohys. Lett. 58 182 
[9] Goldhirsch I and Zanctti G, Clustering instability in dissipative gases, 1993 Phys. Rev. Lett. 70 
1619 

[10] Sela N and Goldhirsch I, Hydrodynamic eguations for rapid flows of smooth inelastic spheres, to 
Burnett order, 1998 J. Fluid Mech. 361 41 

[II] Baldassarri A, Marconi U M B, and Puglisi A, Cooling of a lattice granular fluid as an ordering 

process, 2002 Phys. Rev. E 65 051301 
[12] van Noijc T P C, Ernst M H, Brito R, Orza J A G, Mesoscopic Theory of Granular Fluids, 1997 
Phys. Rev. Lett. 79 411 

[13] Brey J J, Moreno F and Ruiz-Montcro M J, Spatial correlations in dilute granular flows: A kinetic 

model study, 1998 Phys.Fluids 10 2965 
[14] Soto R, Piasecki J and Mareschal M, Precollisional velocity correlations in a hard-disk fluid with 

dissipative collisions, 2001 Phys. Rev. E 64 031306 
[15] Brey J J, Garcia de Soria M I, Maynar P and Ruiz-Montero M J, Energy fluctuations in the 

homogeneous cooling state of granular gases, 2004 Phys. Rev. E 70 011302 
[16] Bobylev A V, Exact solutions of the Boltzmann eguation, 1976 Sov. Phys. Dokl. 20 820 
[17] Pagnani R, Marconi U M B and Puglisi A, Driven low density granular mixtures, 2002 Phys. Rev. 

E 66 051304 

[18] Marconi U M B and Puglisi A, Steady-state properties of a mean-field model of driven inelastic 

mixtures, 2002 Phys. Rev. E 66 011301 
[19] Ernst M H and Brito R, Driven inelastic Maxwell models with high energy tails, 2002 Phys. Rev. 

E65 040301 (R) 

[20] Santos A and Ernst M H, Exact steady-state solution of the Boltzmann eguation: A driven one- 
dimensional inelastic Maxwell gas, 2003 Phys. Rev. E 68 011305 

[21] Cecconi F, Diotallcvi F, Marconi U M B and Puglisi A, Fluid-like behavior of a one- dimensional 
granular gas, 2004 J. Chem. Phys 120 35 

[22] Visco P, Puglisi A, Barrat A, van Wijland F and Trizac E, Energy fluctuations in vibrated and 
driven granular gases, 2006 Eur. Phys. J B 51 377 

[23] Brey J J, Dufty J W, Ruiz-Montero M J, Linearized Boltzmann Eguation and Hydrodynamics for 
Granular Gases, (2004) in Granular Gas Dynamics - Lecture Notes in Physics, Spinger (Berlin) 



